Modeling the quantum evolution of the universe through classical matter 
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We show that the quantum behavior of the radiation-dominated flat Friedmann- Robertson-Walker 
(FRW) universe can be reproduced by classical cosmology given that the universe is filled with an 
exotic matter. In the case of a fiuid, we find an implicit equation of state (EoS) which asymptotically 
approaches p = p/3, showing consistency with the classical model. In the case of two non-interacting 
scalar fields, where one of them is of the phantom type, we find their potential energy. These scalar 
fields behaves asymptotically as radiation. 
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I. INTRODUCTION 

Classical singularities are always present in cosmologi- 
cal models [1] . This is a burden of general relativity and 
there is a general agreement that a full quantum theory 
of gravitation will solve this problem, cither indicating 
how to deal with the singularities or excluding them at 
all. There are several pieces of evidences suggesting that 
quantum mechanics will be able to exclude the classical 
singularities on cosmological models. It was shown in 
[2| 3J that the quantum evolution of the FRW universe 
filled with a perfect fluid with EoS p — wp is non-singular 
and matches the classical evolution when t — >■ oo. In this 
context we have shown that the energy conditions for 
matter are quantum mechanically violated [4], a neces- 
sary requirement for the exclusion of singularities. 

Here we consider another possibility, a classical exotic 
matter which, through Einstein equations, reproduces 
the quantum behavior of the universe. For simplicity, we 
will consider the quantum flat FRW universe filled with 
radiation. We then show that in the classical limit, the 
matter content which fills the universe tends to become 
radiation with EoS p = p/3. This indicates that it is pos- 
sible to get important results from quantum cosmology in 
the classical scenario and still recover the main features 
of classical cosmology as time flows. We will consider two 
kinds of matter. A perfect fluid, for which we will find an 
implicit EoS and two non-interacting scalar fields, one of 
them of the phantom type, for which we will be looking 
for their potential energy. 

Non-linear EoS have already been considered in previ- 
ous works. For instance, in [SI [B] the general quadratic 
EoS p = po + ap -f Pp^ was studied and it was shown that 
with the right choice of the parameters po j ct and /3 this 
exotic fluid can model dark matter. In [B], it was shown 
that such EoS fits a realistic model where a decelerating 
universe turns into an accelerating one. More general 
EoS were studied in [BHS] and again it was demonstrated 
that cosmic acceleration can be modeled by exotic fiu- 



ids. In particular, in [7j, it was proved that the early 
time infiation and late time acceleration can be unified 
by such fluids. An example of implicit EoS can be found 
in [9|. By considering an implicit EoS depending also 
on the Hubble parameter H = a/ a, a universe with one 
big bang and one big rip singularity was found and it was 
shown that the phantom divide w = — 1 is crossed thanks 
to the effect of the inhomogeneous term H in the EoS. 

A universe filled with multiple (canonical and/or phan- 
tom) scalar fields was studied in ^10]. There, they re- 
constructed the entire evolution of the universe, includ- 
ing early inflation and late-time acceleration, using these 
fields. It is also possible to reconstruct the evolution of 
the universe through modified theories, for example F{R) 
gravity. This can be seen in |11| . 

In this paper we search for an exotic matter in a non- 
singular universe. Since the energy conditions are vio- 
lated due to the exclusion of singularities, cosmic accel- 
eration comes as a bonus. However, the expression for 
the scale factor is algebraic in time, so this cosmic accel- 
eration is not enough to generate inflation. 



II. QUANTUM REPRESENTATIVE OF THE 
SCALE FACTOR 



The action for general relativity with a perfect fluid 
in Schutz formalism |12l I13j is given by (in units where 
IQttG = I) 



Sr 




2 / (fxVhhabK"'' 

dM 



(1) 



where hab is the induced metric over the boundary dM. 
of the four-dimensional manifold M, K"'^ is the second 
fundamental form of the hypersurface (9A^ and p is the 
fiuid pressure, which is linked to the energy density by 
the EoS p = wp. 

The super-Hamiltonian for this action in a fiat FRW 
universe given by the metric 
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ds' 



-N^{t)dt^ + a'^{t){dx^ + dy^ +dz^) (2) 



and filled with radiation, is given by 



III. EXOTIC EQUATION OF STATE 



n2 



0, 



(3) 



where Pa = —12^a/N is the momentum conjugated to 
the scale factor a{t) and pr is the momentum conjugated 
to the dynamical degree of freedom of the fluid 14] . As 
we can see in Eq. (|3|, this super- Hamiltonian is a con- 
straint. Following Dirac algorithm !15] for quantization 
of constrained Hamiltonian systems, we make the substi- 
tutions 



Pa 



.d_ 
da 



Pt 



.d_ 
''df 



(4) 



and demand that the super-Hamiltonian operator anni- 
hilate the wave function of the universe. The result is 
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0, 



(5) 



where t — — T is the time coordinate in the conformal- 
time gauge N{t) = a{t) [2]. This is the so called Wheeler- 
DeWitt equation of the universe [16]. Note that in this 
case we do not have operator ambiguities. 

The scale factor is defined on the half-line (0,oo), 
therefore with the usual scalar product 



(*l$) = 



*(a, i)*$(a,t)(ia. 



(6) 



a boundary condition in a = is necessary in order to 
ensure self-adjointness of Eq. (Is]). Lemos found in [2] 
that all possible boundary conditions are given by 



*'(0,i)=/3vl/(0,t), 13 e 



(7) 



For the sake of simplicity, let us consider only two pos- 
sibilities (3 = and (3 = oo, corresponding to Neumann 
and Dirichlet boundary conditions, respectively. Eq. ([5|) 
can be easily solved by finding its propagator [2] or sepa- 
rating variables [3] . Then the expected value of the scale 
factor can be readily computed from the wave function 
^(a,t) through the formula 



(a) it) = 






The result is H |3] 

a(i)=aov/7^+t2. 



(8) 



(9) 



where oq and 7 are positive constants, 7 being related 
to the initial wave packet. Note that when t ^ 00 we 
have that a{t) ~ t. In the cosmic gauge, t is given by 
a(t)dt = dr. Therefore a^r) ^ t^^^, as expected from 
the classical model |T7]. 



The four-velocity of a perfect fiuid in the comoving 
coordinates in the conformal gauge N{t) — a{t) is given 

by 



f/^ = (a(i), 0,0,0), 



(10) 



in such a way that the energy-momentum tensor of the 
fluid 



is given by 



Tpj. = {p + v)U^iUi, +pg^u 



T^,„ == a'^{t)diag{p,p,p,p). 



(11) 



(12) 



Einstein fleld equations (in units where IGttG — 1) 
become 



G, 



00 



= 3^ = ^«^(t)p 



(13) 



Substituting Eq. ^ into the above equations leads to 
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(14a) 
(14b) 



Prior to flnd the EoS of this fluid, let us study the 
expressions for p and p as functions of time. 

In Fig. fll-a) we see that p{t) does not start with a 
singularity. This is expected since the universe with the 
scale factor ^ is non-singular. Then p{t) reaches a max- 
imum pmax = 8/(9aQ7^) at t = 'y/V2, thereafter it starts 
to decay to as i -^ 00 (as the universe expands). Note 
that p(t) is not a one-to-one function, so if we wish to in- 
vert i as a function of p, we must do this in two branches 
(more on this later). We know that the energy con- 
ditions, along with some reasonable restrictions on the 
causal structure of the spacetime, lead to singularities in 
cosmological models. Therefore, for a non-singular uni- 
verse, we expect the violation of the energy conditions. 
This can be verified in Fig. fUb), where we analyzed the 
strong energy condition, which states that p+3p > 0. For 
small t, we have p + 3p < 0, but we recover the classical 
inequality as time flows. Finally, we see in Fig. Ill-c) that 
the classical expression p = p/3 is recovered as t — >■ 00. 
This is not a surprise since a{t) ^ t in this limit. 

We can now find the exotic EoS of the fluid which will 
reproduce the quantum behavior of the universe. By Eqs. 
( 14a I and ( 14b ) we have 
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-P] - iY- 
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(15) 
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Eq. ( 14b ) asserts that 
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FIG. 1; a) Letting ao = 7 = 1, we see that the curve of p{t) 
starts at 0, it has a maximum Pmax ~ 9/8 at i = l/\/2 and 
it decays to as i — >■ cxa. 6) Also, looking to the graphic of 
p{t) + 3p{t) for the same values of ao and 7, we clearly see 
that the strong energy condition is violated for small values 
of t. c). Finally, we consider the graphic of p(t)/ pit) for 
ao = 7 = 1. In the classical limit (t — >■ 00), p/ p — >■ 1/3. 



Adding both equations leads to 



TiP'P)- ,., 



{l^+t^f 



(16) 
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so that 



iip-p) 



4/3 
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2/3 
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(17) 



(18) 



Therefore, we have the implicit EoS f{p,p) = 0, where 



/(P,p) = (f-p) -^(P"P) 



(19) 



Eq. ( |19[ ) can be used to express p as a function of p. 
Doing this inversion we come out with three branches, as 
can be seen in Fig. [2] Two of these branches (denoted by 
dashed lines) are smoothly connected and they represent 
the physical "trajectory" in the diagram p-p. The third 
branch (the dotted hne) is spurious. In fact, in Eq. ( 14a) 
we can invert i^ as a function of p for < p < Pmax- 
We find three different expressions for P (since we are 
working with a third order polynomial in t^), two posi- 
tives and one negative. The negative one is non-physical 
since it corresponds to an imaginary time. Because p{t) 
is not a one-to-one function, it must be inverted in two 
steps. Therefore, one of the positive expressions for t^ 
corresponds to the phase where p grows up to Pmax while 
the other one corresponds to the phase where p decays 
from yOmax to 0. Finally, replacing the expressions of i^ 
in Eq. (14b) we find the three branches of p(p), the one 



corresponding to the imaginary time being spurious. 
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FIG. 2: The three branches of p = p{p) obtained from Eq. 
(191. Two of them (denoted by dashed lines) are connected 
smoothly and represent the physical trajectory of p as a func- 
tion of p as time flows (see the arrows representing the direc- 
tion of time). The other one is spurious. Inset we see that as 
i — >■ 00, we recover the classical EoS p = p/3. 

Now we will show that the classical EoS p = p/3 is 
recovered when t — >■ (X), or similarly, p — >■ and p — >■ 
(see Fig. [2]). To do this, we expand f{p,p) in Taylor 
series up to second order. The result is 



/(p,p) = 2(^-p)Vo(3). 



(20) 



When p — ^ and p ^ we have f{p,p) = ^ p = 
p/3. This shows that when i — >■ oo the classical limit is 
recovered. 



it is easy to show that [by using Eq. ( 24b I 



V{(j)st) + C/(0ph) = -^-^ exp 



-V2{(l)st - </'Ost) ■ (27) 



IV. SCALAR FIELDS 



Equivalently, we have 



We can try to reproduce the quantum behavior of the 
universe by adding a scalar field instead of a perfect 
fluid. We choose to add a phantom and a standard scalar 
field. As we will see in a moment, with these two non- 
interacting fields we can fit the initial cosmic acceleration 
and still recover the classical behavior of the universe as 
t — > oo. 

The Lagrangian for these two fields is given by 



+ ^3'^"a^'/'ph9.</'ph - C/(0ph), 



(21) 



where 4>st corresponds to the standard scalar field while 
(/)ph corresponds to the phantom one. 

The energy density and the pressure for these fields 
can be expressed as 



i?t 



bl^ 



P 



2a? 2a2 
2a2 



2a2 



+ F(0st) + t/(0ph) 
- F(0st) - t/(0ph). 



(22) 



By adding and subtracting the above equations we 
have 



'V = 



i?t 



-a ■ 



p-p = 2T/(0st) + 2C/(</.ph). 



(23) 



Using Eqs. (14a) and (14b) leads to 

2^2- 



v{<i>. 



t - V^ph 

C/(</'ph) 



= 4 



7 



(7^+^2)2 

2 

a2(72+i2)2- 



(24a) 
(24b) 



We demand the standard scalar field and the phantom 
one to take care of the positive and negative parts of the 
right side of Eq. ( 24 ) , respectively. In this way 



t = 7tan('^'^'^), 



so that 

F(0,t) + f/(0ph) 



^h' 



1 + tan2 



/ 0ph-0Oph \ 

V 2 ; 



(28) 



(29) 



Obviously the right hand sides of Eqs. (27) and (29) 



are equal. Therefore, both equations are satisfied if 



V^(4t) = « 



^-^exp 



C/((/.ph) = (l-a)- 



-\/2((/)st - (/-Ost) 

2 



(30) 



<l' 



1 + tan' 



/ 0ph 



— f/'Oph 



where < Q^ < 1- 

It is also simple to check that in this case p — pjZ as 

t — )■ oo. 



V. FINAL REMARKS 

We have found an exotic classical fluid which repro- 
duces the quantum behavior of the universe. This fluid 
obeys an implicit equation of state. When we use this 
equation to find p as a function of p we find three dif- 
ferent branches. Two of them correspond to the physical 
"trajectory" in the p-p diagram. The third one is spuri- 
ous. Moreover, this equation of state approaches p = pjZ 
as t — >■ oo (or p and p goes to zero as the universe ex- 
pands). We can also reproduce the quantum behavior 
of the universe with the insertion of two non-interacting 
scalar fields, one of them being usual and the other of 
the phantom type. We are free to choose the fraction of 
potential energy for each field. In particular, we can set 
one of them to be free. In either case, we showed that 
it is possible to get important results from quantum cos- 
mology and still recover the classical limit with the use 
of classical matter. 
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Since 



7^ -I- 1^ = 7^ exp 



0st(O='^Ost + \/21n(l + tV7') 

</'ph(t) = 0oph + 2 arctan (i/7)- 

(25) 

(26) 
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